Low-energy effective theory of the toric code model in a paraUel field 
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We determine analytically the phase diagram of the toric code model in a parallel magnetic field 
which displays three distinct regions. Our study relies on two high-order perturbative expansions 
in the strong- and weak-field limit, as well as a large-spin analysis. Calculations in the topological 
phase establish a quasiparticle picture for the anyonic excitations. We obtain two second-order 
transition lines that merge with a first-order line giving rise to a multicritical point as recently 
suggested by numerical simulations. We compute the values of the corresponding critical fields and 
exponents that drive the closure of the gap. We also give the one-particle dispersions of the anyonic 
quasiparticles inside the topological phase. 
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I. INTRODUCTION 

For more than 30 years, lattice gauge theories have 
been the subject of intense researches especially in high- 
energy physics where they aim at describing quark 
confinement^. Such theories are deeply related to topo- 
logical phase transitions characterized by the absence of 
local order parameters^'' . One of the most famous mod- 
els where such transitions occur is the Z2 gauge and mat- 
ter theory whose phase diagram in the three-dimensional 
case (2-1-1) has been widely studied by means of various 
methods'''''*'''''*'''. Interestingly and very recently, this 
model has been shown to be equivalent to the toric code 
model (TCM) in a magnetic field in which one introduces 
ancillary (matter) fields together with a gaugc-invariance 
constraint'". 

The TCM was introduced by Kitaev to perform topo- 
logical quantum computation". This spin model can be 
solved exactly and exhibits two kinds of dispcrsionless 
excitations, called charges and fluxes, which have mutual 
anyonic statistics although each of them are bosons. In 
the absence of an external magnetic field, these anyons 
are localized on the vertices (charges) and on the pla- 
quettes (fluxes) of a square lattice. Let us emphasize 
that the detection of anyonic statistics in the TCM has 
been the subject of several experimental proposals in 
optical lattices'^''^'*, although there the TCM appears 
as a low-energy effective theory of Kitaev's honeycomb 
model"'' '. 

The aim of this paper is to study the influence of a mag- 
netic field in the TCM. Contrary to a recent study'", we 
directly consider the quantum problem instead of using 
its classical counterpart. As we will see, the magnetic 
field gives rise to a nontrivial phase diagram which dis- 
plays first-order and second-order transition lines merg- 
ing in a topological quantum multicritical point located 
at the confluent of topological and ordered phases. Ad- 
ditionally, we provide a quasiparticle (QP) description of 
the anyonic excitations in the topological phase. 

To investigate this issue, we use several perturbative 



treatments. First, we perform a standard (linear) spin- 
wave analysis'^'' which captures quantum fluctuations 
around the classical ground state and is thus certainly 
valid (qualitatively) in the large- field limit. Second, we 
compute the perturbative expansion of the ground-state 
energy and the gap in the small-field limit and in the 
large-field limit by means of the continuous unitary trans- 
formations method' ''^*''^'''^". These three approaches al- 
low us to propose a consistent picture of the phase dia- 
gram. Finally, we also give the dispersion relation in the 
1-QP sector in the topological phase. 



II. MODEL 

We consider the following Hamiltonian 
i7=-J,^^-Jp^Sp-/^,^<~;^,^a^ (1) 

s p i i 

where the af's are the Pauli matrices. As = Yiies'^f^ 
and Bp = Yii^p '^i ■ Subscripts s and p refer respectively 
to sites and plaquettes of a square lattice whereas i runs 
over all bonds where spin degrees of freedom are located 
(see Fig. 1). 

Up to a global normalization, the parameter space of 
Hamiltonian (1) is three-dimensional. Here, we focus on 
the two-dimensional subspace defined by Js = Jp = J 
which, for h^ — h^ = 0, coincides with the TCM". We 
emphasize that this subspace is not the same as the one 
considered in Ref.'" where the variables {Js, Jp,hx,hz) 
are linked via the mapping onto the isotropic Z2 gauge 
Higgs model. Consequently, one cannot compare our re- 
sults with the numerical data'" but, as we shall sec, our 
phase diagram displays similar qualitative features. 

Let us mention that the single-component magnetic 
field case has also been addressed recently'^''^'^ but its 
low-energy properties are exactly the same as those of 
the celebrated two-dimensional Ising model in a trans- 
verse field whose phase diagram has been determined ac- 
curately many years ago^\ 
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FIG. 1: (Color online). A piece of the square lattice on which 
the TCM is defined. Spins (dots) are located on the bond and 
interact if they share either a vertex (s) or a plaquette (p). 



III. LINEAR SPIN- WAVE THEORY 

As a first approach, let us compute the low-energy 
spectrum of H in the semi-classical (large-spin) limit. 
Therefore, we consider the Hamiltonian Hs obtained by 
replacing the Pauli matrices by SU(2) spin-S operators 
(af -^ Sf/S) which are considered as classical vectors in 
the large-S* limit [Si = 5'(sin6'i cos(/)i,sin0i sin(/),i, cos6'j)] 
. Assuming that the ground state is given by a uniform 
configuration {0i = 9,4>i = 0, Vi)~"\ we are led to min- 
imize the corresponding classical energy E. Next, we 
use the Holstein-Primakoff representation of the spin op- 
erators to compute the leading 1/ S corrections and we 
obtain the excitation energies A^ = i -^/e^ — 7^, where 

Efe = 7fc-f 2j(sin^ 6'0 + cos'* do^+hx sm 9a + hz cos^O: (2) 

and 

7fc = — Jeos^ 6*0 sin^ ^0*^ 2cos(fe.ni) -I- 2cos(fe.n2) 

-I- cos [fc.(ni + 722)] + cos [k.{ni - 722)] >. (3) 

Here, rii and n2 are unit vectors displayed in Fig. 1 and 
9o is the angle which minimizes E. 

Within this approach we find that the low-energy spec- 
trum is gapped for all values of the field except for 
hx = hz = y/2J = h™ where it vanishes. Further, we 
find that below this value, on the isotropic line h^. = hz, 
the classical energy has two degenerate minima with the 
same characteristics. This leads us to conclude that the 
segment < h^ = h^ < h^™ is a first-order transition line 
which ends at the second-order critical point h^™ . Note 
that the value of h™ , being independent of S, certainly 
differs from the actual value for S = 1/2. 

Although one can reasonably believe in the qualitative 
features of the spin- wave scenario in the large-field limit, 
it certainly fails for small fields. Indeed, in this limit, 
the ground state is far from being a separable state for 
5 — 1/2 as can be inferred from the zero-field case. In 
addition, for S > 1/2, Hs does not commute with A's 



and B's in zero field which is a key ingredient of the 
TCM's topological character. 



IV. LARGE-FIELD LIMIT {h^,h^ > J) 

To determine the value he of the critical field for 
S — 1/2, let us investigate the low-energy spectrum 
of Hamiltonian (1) in the strong- field limit and on the 
isotropic line h^ = hz = h. For J = 0, the ground state 
of H is fully polarized in the field direction, and elemen- 
tary excitations are static single spin-flips with energy 
cost 2^^^h. For J > and setting h = 2~^/^, Hamilto- 
nian (1) can be recast into 



H = -N 



E 

n=0,±l,±2,±3,±4 



Tn 



(4) 



where the operators T„ are proportional to J and change 
the number of excitations Q by n, i.e., [QjTn] = nTn- 
Their expressions are easily obtained but are too lengthy 
to be given here. To study this Hamiltonian, we 
used the perturbative continuous unitary transforma- 
tions (PCUT) method^''*'^"'-", which allows one to con- 
struct, order by order (in J), an effective Hamiltonian 
which is unitarily equivalent to H but conserves the num- 
ber of QPs. The ground state of the effective Hamiltonian 
is the 0-QP state, whereas the lowest-excited states lie in 
the 1-QP sector. We thus have access to the ground- 
state energy per spin cq, to the dispersion of the QP and 
consequently, to the gap A. At order 5 in J, one has 
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A standard Dlog Pade approximants analysis for the 
gap leads to he = 0.48(2)J which, as anticipated, 
strongly differs from the spin-wave value {h™ — \/2J) 
but confirms the existence of a critical point on the 
isotropic line and undoubtedly provides a more accurate 
determination of its location. Moreover, in the vicinity of 
this point and on this line, a [2, 2]-Dlog Pade approximant 
yields A ~ (ft, — h^Y with v ~ 0.73 instead of z^ = 1/2 
in the linear spin-wave approach. This result is hardly 
compatible with an Ising-type critical point suggested in 
the 2,2 gauge Higgs model^ for which v = 0.6301(8)-"', 
but this is clearly due to the relatively low order of our 
expansion. 



V. SMALL-FIELD LIMIT {h^,h^ < J) 

This region is, by far, the most interesting and chal- 
lenging one. As explained above, for h^ = hz = 0, Hamil- 




FIG. 2: (Color online). One-flux (up) and one-charge (down) 
dispersions for h^ = 0.1, hx ~ 0.05, and J — 1/2 in the 
reciprocal lattice. In this case {hx < h^), the gap is given by 
the minimum of the one-charge dispersion [A ~ 0.54 as can 
be checked from Eq. (9)]. 



tonian (1) is exactly the TCM and one thus expects the 
system to be in a topological phase, whose breakdown 
should occur at finite fields only. The ground state of 
the TCM has all Ag = Bp = 1. Elementary excitations 
are such that one charge or one flux is present, namely 
one As = —1 or one Bp = — 1 , and have an energy cost of 
2 J. Note that here, we are assuming open boundary con- 
ditions and the thermodynamic limit which allows one to 
consider a single (flux or charge) excitation. On a torus, 
such conflgurations are prohibited and one must analyze 
states with pairs of fluxes or of charges. 

In the TCM, fluxes and charges are static but, for fi- 
nite fields, these excitations acquire some dynamics and 
can be considered as true QP. Of course, setting up such 
a description in this (liquid) topological phase is a highly 
nontrivial task. In the present case, this is made possible 
thanks to the strict locality of the anyonic excitations for 
hx ~ hz = 0. Thereafter, we establish perturbatively 
this QP picture and we present results for the one-flux 
and one-charge dispersions in the topological phase. The 
study of the one-particle gaps will allow us to determine 
the boundaries of this phase and to compute the corre- 
sponding critical exponents. 

For nonvanishing fields and setting J = 1/2, the 
Hamiltonian can be written as 



where now Q counts the total number of charges and 
fluxes. The T„ operators are linear in hx and hz and 
satisfy as previously [Q,T„] = nT„. Their precise ex- 
pressions do not bring any special physical insight and 
are thus omitted here. 

Once again, such a form is well suited to a PCUT treat- 
ment. We emphasize that our study amounts to comput- 
ing transition amplitudes of the effective QP-conserving 
Hamiltonian between the highly entangled eigenstates 
of the TCM^^. To this end, it is essential to keep 
track of the anyon positions and of the underlying spin 
background simultaneously. This makes our perturba- 
tion theory more complicated than the one derived in 
Rcf. whose unperturbed Hamiltonian corresponds (in 
the gauge-theoretical reformulation used in Ref.^"), to 
Jp = hz = 0, which has separable eigenstates. 

We have obtained the ground-state energy per spin 
cq, as well as the dispersions of the QP (dressed charges 
and fluxes) which are obviously mapped one onto the 
other when exchanging h^ and h^- A typical dispersion 
is displayed in Fig. 2. The gap A is the minimum of both 
dispersions, and we give it here in the region hx ^ hz 
where charges are the lowest-energy excitations. Both 
eo and A were computed at order 8 in {hx,hz), and we 
obtained 
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H 



-N + Q + To + T+2 + T^2, 



(7) 



When one of the magnetic field components vanishes, 
the TCM is equivalent to the Ising model in a transverse 
field so that setting hx = 0, we recover the results ob- 
tained by He et al?^ . Equation (9) predicts a gap which 
vanishes continuously along a line starting from the crit- 
ical point [hx = 0,hz = 0.1642(2)] and ending at the 
multicritical point [hx = hz = 0.1703(2)]. As previously, 
error bars are obtained from various Dlog Padc approx- 
imants used to analyze the series. We also computed 
the critical exponent driving the closure of the gap. We 
found Vh^=o — 0.65, close to the expected Ising expo- 
nent and Vh-,=h, — 0.70 at the multicritical point (note 
that at order 8, these exponents are not yet fully con- 
verged). For intermediate values, the exponent sticks to 
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FIG. 3: (Color online). Phase diagram in the plane {hx,hz) 
for J = 1/2, where second-order transition lines are drawn 
full (red), the first-order transition line is dashed (blue), the 
multicritical point is a full circle and the critical point is rep- 
resented as an empty circle. Phases I, II and III are discussed 
in the text. 



the Ising value except in the vicinity of the multicritical 
point, indicating that phase transitions arc Ising-likc, ex- 
cept at the multicritical point. This is confirmed by our 
QP picture. Away from this point, only one kind of par- 
ticle condenses (charge or flux), the other one remaining 
gapped. At the multicritical point, both types of parti- 
cles condense simultaneously, and their mutual semionic 
statistics should become important and gives the transi- 
tion an unconventional character. 



VI. DISCUSSION 



where (Ag) ~ 1 and (Bp) ~ 1 in the ground state 
for hx,hz <C J. This phase has dispersive charge and 
flux excitations. Phase II (III) is such that {af) ~ 1 
if /ij; > /iz > J ((ct|) ~ 1 if /i^ > /ij, > J), and these 
phases have dispersive spin-flip excitations. However, us- 
ing Eq. (5) and the Hellmann-Fcynman theorem for the 
ground state energy one can compute (erf) and (erf) in 
phase I, and check that they do not vanish. 

Actually, no local order parameter can be used to char- 
acterize these various phases showing that the previous 
description is very rough. Furthermore, simple non-local 
order parameters can only be found on the Ising lines 
{e.g. if hx = Q, & semi-infinite string of cr^ operators 
on a line of the square lattice of Fig. 1). Finding order 
parameters for these phases thus remains challenging. 

Apart from the phase diagram, the central result of our 
work is the set up of a QP picture for fluxes and charges 
in the topological phase. As a consequence, we have been 
able to compute the phase boundaries and the critical ex- 
ponents by studying the locus of points where the 1-QP 
gap vanishes. This QP description offers a wide range of 
perspectives. Indeed, the present PCUT approach is par- 
ticularly well adapted to study the many-QP physics"''. 
This will allow us to investigate the likely existence of 
bound states in this model made up of flux-charge com- 
posites (fermionic statistics) that may change the critical 
properties. Such fermions are expected to play a major 
role when switching on a magnetic field in the y-direction 
since at lowest order, such a field induces a hopping of 
these fermions (as well as a local transmutation of two 
charges into two fluxes and vice- versa) . Note that it may 
be difficult to study this problem using Monte Carlo sim- 
ulations because of the usual sign problem, which arises 
in the presence of a transverse field. 



The phase diagram obtained from our analytical cal- 
culations is shown in Fig. 3. Second order transition 
lines are obtained from the small-field expansion, the 
first-order line from the classical (large-5) analysis and 
the position of the critical point from the large-field ex- 
pansion. As in Ref.^", phase I is a topological phase 
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